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Abstract 

We study the homogenization of an obstacle problem in a perfo- 
rated domain, when the holes have random shape and size. The main 
assumption concerns the capacity of the holes which is assumed to be 
stationary ergodic. 



1 Introduction 

Let (Q, T, V) be a given probability space. For every wel! and every e > 0, 
we consider a domain D e (uj) obtained by perforating holes from a bounded 
domain D of M. n . We are interested in the asymptotic behavior as e — ► of 
the solution of the following obstacle problem: 

min | J ^|Vii| 2 - fudx; u>0 a.e. in D \ D £ , uE Hq(D) 

for some / £ L 2 (D). This is a well known homogenization problem and the 
asymptotic behavior of the solutions strongly depends on the size and the 
repartition of the holes 

T e {u) = D\D £ . 

This problem was first studied in the case of periodic domains by L. Carbone 
and F. Colombini [CC80] and then in a more general framework by E. De 
Giorgi, G. Dal Maso and P. Longo [DGDML80| and G. Dal Maso and P. 
Longo [DML81], G. Dal Maso [DM81]. Our main reference for this work will 
be the papers of D. Cioranescu and F. Murat |CM82a[ ICM82b] . in which the 
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case of a periodic repartition of the holes D\D £ is studied. It is proved that 
when the number of holes and their size are evolving in a critical fashion, 
then the limiting problem is no longer an obstacle problem, but a simple 
elliptic boundary value problem with a new term that takes into account 
the effect of the holes. 

Our goal is to generalize their result to the case where the holes are 
still located in small neighborhoods of the points of the lattice eZ n but 
have random size and shape. More precisely, we assume that for any e and 
u> the domain D £ (uj) is obtained from a fixed set D by perforating holes 
{S £ (k,uj); k G Z n } such that 

S £ (k,u) C B £/2 (sk) forallfceZ™. 

We denote by 

T £ (iv) = U ke znS £ (k,Lj) DD 
the union of all the holes in D. We then have 

D £ (u) = D(u)\T £ (u). 

The assumptions on the sets S e (k, u>) will be made precise in the next section. 
We can already point out the fact that we will not exclude the case where 
S £ (k,uj) = for some k, thus allowing the fact that no holes may be present 
at some lattice points. 

With these notations, we rewrite the obstacle problem as follows: 

f{ut)= inf J{v\ u £ eK £ (1) 

v£K e 

with 

</(«) = f \\Vu\ 2 -fudx 

JD 1 

and 

K £ = {v G Hl(D) ;v>0 a.e. in T £ }. 

Since K £ is closed, convex and not empty, ([1]) has a unique solution u £ G K £ . 
Moreover, u £ solves 

' -Au e = f in D £ 

< u £ (x) > on T £ (2) 
k u e {x) = on dD\T £ 
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As mentioned in the introduction, it is expected that under appropriate 
assumptions on the size of the holes S e (k,u), the function u e converges 
weakly in H 1 to u solution of 



where U-(x) = max(0, —u(x)). 

The assumptions and the result are made precise in the next section. 
The proof of the main theorem, which is details in Section [3j relies on the 
construction of an appropriate corrector. This construction is detailed in 
Sections U] and [5j first in the case where the holes are balls in dimension 
n > 2, then when no assumptions are made on the shape of the holes (in 
dimension n > 3 only). 

2 Assumptions and Main result 

First, we need to make precise our assumptions on the holes S £ (k,uj). The 
first assumption is mainly technical: 

Assumption 1: There exists a (large) constant M such that for all 
k € Z7 1 and a.e. oj € Q we have 

S e (k, uj) C B M£ n/(n-2) (ek) if n > 3 
S £ (k, uj) C -B e xp(- M £-2) (ek) if n = 2 

for e small. 

As mentioned in the introduction, the asymptotic behavior of the u e 
strongly depends on the size of the holes. The critical size for which inter- 
esting phenomena is observed corresponds to finite, non trivial capacity of 
the set T £ . More precisely, we assume: 

Assumption 2: For all k G IT 1 and a.e. uj € there exists j(k,oj) 
(independent of e) such that 



Au — oiqu 



f mD 



u = 



on dD. 



cap(S £ (k,uj)) =e n j(k,io) 



where cap (A) denote the capacity of subset A of M n , defined by: 



cap (A) = inf 



{ 




V/i| 2 dx ; h G H l (R n ), h > 1 in A, lim h(x) = 



3 



in dimension n > 3 and by 



cap (^4) = inf 



\[ \Vh\ 2 dx; h G H%(Bi), h > 1 in A 



UBi 



in dimension n = 2 and for sets A C B\. Moreover, we assume that there 
exists a constant 7 > 0: 



Finally, our last assumption will be necessary to ensure that some aver- 
aging process occur as e goes to zero: 

Assumption 3: The process 7 : 7L n x f2 1— > [0, 00) is stationary ergodic: 
There exists a family of measure-preserving transformations : 0, — > f2 
satisfying 

7(/c + A/, cj) = 7(/c, Tyoj) for all fc, A;' G Z n and w G fi, 

and such that if i C and r fe A = ^4 for all G IT 1 , then P(A) = or 
P(A) = 1 (the only invariant set of positive measure is the whole set). 

Let us make a few remarks concerning those assumptions: First of all, we 
stress out the fact that the shape of the holes S e is left unspecified and may 
change with e; Only the rescaled capacity is independent on e. The first 
assumption, which implies that the diameters of the holes decrease faster 
than e, guarantees that the capacities of neighboring sets separate at the 
limit (i.e. that cap(US' e ) ~ ^cap(S' e )). And the choice of scaling for the 
capacity guarantee that cap(T e ) remains bounded as e goes to zero (since 
#{Z ra n e~ 1 D} < Ce~ n ). Finally, the hypothesis of stationarity is the most 
general extension of the notions of periodicity and almost periodicity for a 
function to have some self-averaging behavior. 

Under those assumptions, we prove the following result: 

Theorem 2.1 Assume that n > 3 or n = 2 and the holes are all balls. 
Then there exists Qo > such that when e goes to zero, u £ converges weakly 
in H 1 to a function u solution of the following minimization problem 



j(k, to) < 7 for all fceZ" and a.e. wefi. 



(3) 



mm 



{ 
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where u-(x) = max(0, —u(x)). In particular, u solves 



Am — ocqu. 



f in D 



u = 



on dD. 



Moreover, if there exists 7 > such that 



j(k,u) > 7 for all k € 7L n and a.e. w6D 



then «o > 0. 

The general result holds also in dimension n = 2 when the holes have 
random shape. However, because the fundamental solution of Laplace's 
equation is different in that case, the proof is slightly different and more 
technical. 

As in Cioranescu - Murat [CM82a, CM82bJ, the proof of this result relies 
on the construction of an appropriate corrector. More precisely, the key is 
the following result: 

Proposition 2.2 Under the assumptions listed above, there exists a non- 
negative real number «o o,nd a function Wq(x,u) such that 



' Aw 6 = ao 
< w £ (x) = 1 
k w £ (x) = 



in D e (u) 

in T £ (uS) 

on dD\T £ {u) 



for almost all u> & ft, and 







H (D)-weak a.s. uj E f2. 



Note that as in |CM82aj . the equation 



Aw e = ao in D £ (lu) 



can be replaced by the weaker condition: 



For all sequences v £ satisfying: 
( v e = on T e 
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The proof of Proposition 12,21 will occupy most of this paper. It will be 
split in two parts: In Section^ we consider the (simpler) case when the holes 
S e (k,uj) are all balls of random radius. In Section [5l we will use this first 
result to treat the general case (when the holes have unspecified shapes). 

Before turning to this proof, we briefly give, in the next section the proof 
of the main theorem. 



3 Proof of Theorem [2J] 

First of all, standard elliptic estimates give the existence of a function u 
such that 

u £ — ► u H 1 — weak. 

If we introduce the limit energy 

f 1 1 
/aM = j^-\X7v\ 2 + -a v 2 _- fvdx, 

it is readily seen that all we need to show is the following inequality: 

This relies on the following two lemmas: 
Lemma 3.1 For any ip € Wq'°° , we have 

lim / \Vw £ \ 2 ip alx = / ocopdx 

Lemma 3.2 Ifu £ ^uin H 1 -weak, then 

liminf J {u £ ) > f a {u) 

Let us see that those two lemmas imply the theorem: For any v E 
the function v + V-W £ is non-negative on the holes, and is thus admissible 
for the initial obstacle problem. In particular by definition of u £ , we have 

f{u £ ) < f{v + v-w £ ). 



6 



We write 



v + v-w £ ) = j i[| Vu| 2 + \Vv-\ 2 w e2 + \v-\ 2 \Vw e \ 2 } dx 



+ / \v-X7v-w e Vw e + Vf X7v-w £ + Vvv^X7w £ ] dx 



and it is readily check that Lemma 13.11 and the weak convergence of w e to 
in H 1 implies 

Xm^J(y + v^w £ ) = J^aiv), 
as soon as v € W 2 '°°. We deduce: 

c/a(v) > limsup ^f(u e ) 

for all v E Wq'°° . Together with Lemma 13.21 this gives 

for all v £ Wq' 00 . We deduce Theorem 12.11 by a density argument. □ 

Proof of Lemma l3.lt We recall the proof of Cioranescu-Murat [CM82a] : 
Since 1 — w £ = in T e , we have: 

/ Aw £ f{l-w £ )dx = / (p\Vw £ \ 2 dx- V(pX7w £ {l-w £ )dx 
Jd e Jd e Jd e 

and so 



ao(p(l — w £ )dx = 1 a.Q(p(l — w £ ) dx 
D Jd £ 

= / ip\Vw £ \ 2 dx- Vip\7w £ (l-w £ )dx. 
Jd £ Jd £ 

finally 

/ V(pS7w £ (l - w £ )dx = / S7cpVw £ - VipS/w £ w £ dx — >0 

since w £ goes to zero H 1 -weak and L 2 -strong. The lemma follows. □ 

Proof of Lemma l3.2t See Cioranescu-Murat }CM82bj . Proposition 3.1. 
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4 Proof of Proposition [2T2t Balls of random radius 



Throughout this section, we assume that the sets S e (k,uj) are balls centered 
at ek. Since 



cap(B r 



n(n — 2)u) n r n 2 if n > 3, 
2tt 



log r 



if n = 2 



Assumption 2 becomes in this framework: 

S £ (k,uj) = B a e ir{k)UJ )){ek) for all k G 



with 



and 



a e (r) 



rE »/(n-2) if n > 3, 

exp(— r _1 e -2 ) if n = 2, 



r(k,Lj) = I \n{n-2)u) n ) 

j(k, lv)/2tt 
Note in particular that the process 



if n > 3, 
if n = 2. 



r:Z n xO^[0,oo) 
is stationary ergodic and satisfies 

r(k,co) < r for all G Z n and a.e. w G SI 



(5) 



for some constant r > 0. Without loss of generality, we can always assume 
that r < 1/2 (so that there is no overlapping of the holes for any e < 1): 



4.1 The auxiliary obstacle problem 

After rescaling, we look for the corrector w 6 {x,lo) in the form 

w £ (x,uj) = e 2 v e (x/e,Lo) 

with v £ {y,uj) solution to 



Av = a, 
v = e~ 2 



ins 1 D e , a.e. to G Q 

on Ufc g z™ B^{k^){k) 



with 



and satisfying 



Tf(r) 



■ r£ 2/(n-2) ifn>3, 
e _1 exp(— r _1 e~ 2 ) if n = 2, 

e 2 v £ (x/e) — ► in i/ 1 -weak . 



One of the main tool in the proof is the fundamental solution of the 
Laplace equation, given by: 



1 



1 



h(x) = < 
In particular, we note that 

^ls^(r(fc,u,))(0) = 



n(n — 2)to n \x\ n ~ 
log \x\ 



n{n — 2)uj n r n 2 



if n > 3, 



if n = 2. 



if n > 3, 



2vr 



(log(e) +r- 1 e" 2 ) if n = 2, 



so we expect the rescaled corrector v £ (x, uj) to behave near the hole B^ k u! ^ (k) 
like the function 



h k (x) := < 



r(k, uj 



,n-2 



-y(k,u) h(x-k) = | x ^ fc[ra _ 2 , ifn>3 
7(fc, /i(x — k) = —r(k, uj) log \x — k\ if n = 2, 



where 



j(k,u) 



(r(k,uj)) n - 2 n(n-2)uj n if n > 3 
2irr(k, uj) if n = 2. 

Since /i^ satisfies 

A/ifc = —j(k,oj) 5(x — k), 
we will construct v £ (x,uj) by solving 

J At; = a - Efcez^nA 7(*. v)5(x - k) in D , 
) i> = on 

The main issue is thus to find the critical a for which the solution of the 
above equation has the appropriate behavior near x = k. 
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Following [CSW05] . this will be done by introducing the following ob- 
stacle problem, for every open set A C R n and a € M: 



v a ,A(%, = inf < v{x) ; Av < a — j(k, lu)5(x — k) 

Clearly, the function v a ,A is solution of 

Av = a — j(k,oj)6(x — k) 

kez n nA 

whenever it is positive. Note that the function 



v > in A 
v = on dA 



a 



2n 


\x - 1 








a 




— \x 


< 


2n 




a 








2n ,X 



(6) 
(7) 

(8) 



r(k, uj 



,n-2 



n-2 • 



also satisfies 



\x — k\ 



Ah at k(x) = a — r y(k,uj)5(x — k). 



if n > 3, 



It follows from (JT]) and the maximum principle that if B\(k) C A, then, for 
all x in B\{k) and for almost every u in Q, we have 



a 



h a , k {x) - — - r n - 2 ifn>3 
An 



a 



if n = 2. 



0) 



4.2 Critical a 

The purpose of this section is to prove that for a critical a, v aA behaves like 
h a ^k near S £ (k,u>). For that purpose, we introduce the following quantity, 
which measures the size of the contact set: 

m a (A,uj) = \{x G A; v a:A (x,uj) = 0}| 

where |^4| denotes the Lebesgue measure of a set A. 

The starting point of the proof is the following lemma: 
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Lemma 4.1 The random variable rn a is subadditive, and the process 

T k m(A, uo) = m(k + A, u) 
has the same distribution for all £ Z". 

Proof of Lemma I4.lt Assume that the finite family of sets (Ai)i^j is such 
that 

Aid A for all? £ I 
Ai n Aj = for all i / j 
\A-U ia Ai\ = 

then v a! A is admissible for each Ai, and so Vct,A{ < u a,A- It follows that 
{v a , A = 0} n Ai C = 0} 

and so 

m a (A,u) = ^2 \{vot,A = 0} n Ai\ < ^2 \{V<x,Ai = 0}| = y^fn a (Ai,oj), 
iei iei iei 

which gives the subadditive property. Assumption 3 then yields 

T k m(A,Lo) = m(A,r k u) 
which gives the last assertion of the lemma. □ 

Since m a {A,Lv) < \A\, and thanks to the ergodicity of the transforma- 
tions Tfc, it follows from the subadditive ergodic theorem (see [DMM86 ) that 
for each a, there exists a constant £(a) such that 

m a (B t (o),u) _ - 

\B t (0)\ ~ e{a) a - S '' 

where -Bt(O) denotes the ball centered at the origin with radius t. Note that 
the limit exists and is the same if instead of -Bf(O), we use cubes or balls 
centered at txo for some xq. 

If we scale back and consider the function 

w £ a (y,uj) = e 2 v atB ^ l{£ -i xo) (y/e,uj), in Bi(x ), 

we deduce _ 

,. \{y; w e a (y,u) = o}| -. . 

hm : : = i(a) a.s. 

s-o \Bi\ v ' 

The next lemma summarizes the properties of £(a): 
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Lemma 4.2 

(i) £(a) is a nondecreasing functions of a. 

(ii) If a < 0, then t(a) = 0. Moreover, if the radii r(k,co) are bounded from 
below, then 1(a) = for any a such that a < n(n — 2) inf^g^n r(k, lu) u ~ 2 
almost surely. 

(Hi) If a > 2 n n(n — 2) sup fceZ n r(k, uj) n ~ 2 (or a > 8r for n = 2) almost 
surely, then 1(a) > 0. 

Proof. 

(i) The proof follows immediately from the inequality 

v<y,A < Va',A f° r anv a i a> such that a < a. 

(ii) If a is negative, then the function ^-\x — xo\ 2 — ^(tr) 2 , which is a sub- 
solution of ((7J), is positive in tB r (xo) and vanishes along d(tB r (xo)) for any 
ball B r (xo) and for any t > 0. We deduce: 

Va.iB > - x | ~ ^\ tr ) >0mtB r (x ) 

for all i > 0. Therefore m a (tB,uo) = for all i > 0, so £(a) = for all 
a < 0. 

Furthermore, if r(k,tu) is bounded below: 

r(/c, w) > r > for all k E Z n , a.e. well, 

then, the function ^\x — k\ 2 + npjjn=i ~ ^ ~ r n ~ 2 is a solution of ([7|) in 
B\(k) which vanishes on dB±(k) and is strictly positive in B\(k) as long as 
a < n(n — 2)r n ~ 2 . As above, we deduce that m a (tB,u) = for all t > 
and for all a < n(n — 2)r n ~ 2 . 

(iii) The function h a ^(x) = ^\x — k\ 2 + ^L^n-i is radially symmetric and 
reaches its minimum when 



| x — k\ = R(a, k) := < 



n(n-2)r(k,uj) n - 2 \ 1 / n 



2r(fc,u) 



a 

1/2 



when n > 3 

(10) 

when n = 2 



In particular, for a > 2 n n(n — 2)r(k,u)) n 2 (or n > 8r(k,u> when n = 2), 
we have R(a, k) < 1/2 and so the function 



5fc(a0 



h a>k (x) - D k in B R{aM) {k) 

inR^i^)^) 
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satisfies 

A<?fc < a — j(k, uj)5(x — k) in C\(k) , 

and 

g k = md(k)\B 1/2 (k) 

where C\{k) denotes the cube of size 1 centered at k, and the constant Cfc 
is chosen in such a way that and vanish along dB^i a k y. 

[ | (1 — log(2r/a)) when n = 2. 

By definition of U Qj j_b, we deduce that 

v a ,tB(x) < 9k{x) in i£ a.s. 

In particular, this implies that v at tB vanishes in tB \ Ufc e z n -Bi/2(^)> and so 

m a (tB,uj) f\Ci\-\B 1/2 \\ ljJn 

> rrr-, = 1 — — — a.s. 



\tB\ ~ \ \d\ J 2 

We conclude 



( a ) > i _ > 0. 



□ 



Using Lemma 14.21 we can define 

«o = sup{a ; 2(a) = 0}. 

Note that ao is finite under Assumption 3 (Lemma l4.2l (iii)) and that ao > 
is strictly positive as soon as the r(k,u) are bounded from below almost 
surely by a positive constant (Lemma 14.21 (ii)). 

In the rest of this section, we are going to show that the function 

El \ ■ A I \ A / ■ rn W>l0nT £ flD 

w (x,u)) = mf < w(x); Aw < a m D \T e , «n\^ 

satisfies all the conditions of Proposition 12.21 We will rely on a series of 
intermediate functions. 

For the first lemma, we fix a bounded subset A of W 1 and we denote by 

v £ q (x,uj)=v -i a (x,uj) (12) 
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the solutions of © defined in e 1 A. We also introduce the rescaled function 

w e a {y^) = e 2 v £ a (y/e,uj), 

defined in A. 

The key properties of v £ a are given by the following lemma: 
Lemma 4.3 

(i) For every a and for every fc£Z", we have 



_ r n~2 



K,k{ x ) 



2n 



if n > 3 



v e Jx) > < 



a 




ifn = 2 




if a (x) < h a;k (x) + o(e 2 ) 



for all x € Bi^ik) and almost everywhere to E ft. 



Since 




(r(fc,u)))(0) 



= < 




we deduce the following corollary: 



Corollary 4.4 

(i) For every a and every k € Z n such that r(k,u) > 0, we have 
lf a {x) > e~ 2 + o(l) on dBa^ r (k,u)){k) a.e. wed 



and so 



uf a > 1 + o(e 2 ) on 8T £ (uj) a.e. uj € $7 



for all a. 
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(ii) For every a > cto o,nd every k G Z n ; we have 

V £ a {x) < £~ 2 + o(e~ 2 ) on <9-E>a e (r(fc,u;))(^) a - e - ^£0 

and so 

< 1 + o(l) on dT £ (uj) a.e. uj G O 

Proof of Lemma \4-3[ 

(i) Immediate consequence of ©. 

(ii) 

Preliminary: First of all since ^4 is bounded, we have 

A c Br(z ). 

Without loss of generality, we can always assume that Br(xq) = Bi(0). We 
then introduce 

v%(x,(j) =v a>e -i Bl (x,uj), 

the solutions of © in B e -i(0). It is readily seen that v £ a is admissible for 
([6]) and thus 

v £ a (x,u)) < v e a (x,uj) for all x G a.e. G fi. 

It is thus enough to prove (ii) for v e a . 

We will need the following consequence of Lemma 14.11 (see |CSW05j for 
the proof): 

Lemma 4.5 For any ball B t {xq) G £?i(0), i/te following limit holds, a.s. in 
ijj 

lim |R(x,u;) = 0}n^ £ -i r (g- 1 x 1 )| = - 



Step 1: We can now start the proof: For any 5 > 0, we can cover B £ -i by 
a finite number iV (< C5~ n ) of balls -Bj with radius 5e~ l and center £ —1 ;cj. 
Since a > ctQ, we have Z(a) > 0. By Lemma H3J we deduce that for every 
i, there exists Si such that if e < e%, then 

|{«^(x,a;) = 0} n Bi\ > a.s. u. 

In particular, if e < inf £j, then v e a {yi) = for some yi in Bi a.s. wefi. We 
now have to show that this implies that v e a remains small in each Bi as long 
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as we stay away from the lattice points k G Z n . More precisely, we want to 
show that 

sup v £ a < C5 2 e' 2 . 

Step 2: Let 77 be a nonnegative function such that < r/(x) < 1 for all x, 
rf(x) = 1 in Bi/ 8 and 77 = in R n \ By^. Then the function u = v e a * 77 is 
nonnegative on 2£?j and satisfies 

-C < An < C 

where C is a universal constant depending only on n and f. In particular, 
since Bi has radius Harnack inequality yields: 

supu < Cinf u + Caffe" 1 ) 2 . 



Step 3: We need the following lemma: 
Lemma 4.6 If Av < a in B r (yo), then 



— [ v(x) dx < v(yo) + aC(n)r 2 

B r J B r (yo) 



where C(n) is a universal constant. 

Proof: We note that the function v(x) — j^\x — yo\ 2 is super-harmonic in 
B r (yo). The lemma follows from the mean value formula. □ 

Now, we recall that v^ivi) = and Av% < a in B 1 u(yi). So 

1 



B1/4 JB 1/4 (y x ) 



v e a {x) dx < v £ a {yi) + aC(n) 



In particular, we have 



i(yi) < / v £ a (x)dx <C(a,n) 



Step 4-' Steps 2 and 3 yield 



supu < C(a,n)(l + a(5e^ L Y). 
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and since Av% > in Bi \ n k £Z n {k}, we have: 

1 



B, 



1/8 JB 1/a (y) 



v e a {x) dx < Cu{y) 



for all y £ Bi\ n fe6Z »S 1/4 (fc). 

It follows that for every 5 and for e small enough, we have: 

sup v £ a < C5 2 e- 2 . 

B £ -i\U fcez «S 1/4 (fc) 

The definition of v £ a and the fact that h a }~ > on dB-yji implies that 



v%{x) < K )k (x) + C5 z e 



2„-2 
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By 2 (k) 



for all k 6 Z n . □ 



We now want to use the solution (|12p of the obstacle problem ([6|) with 
A = D to study the properties of the free solution Wq of 

f Aw = ao- EfceZ"nz? 7(*i w)<J(z - efc )> in ^ 
\ = on dD 

We prove: 

Lemma 4.7 For every k € Z n , satisfies 

h e ak {x)-o{l) < w £ (x) < h £ ak (x)+o(l) Vx € B E/2 (ek)<lD a.e. wGO, 

(13) 

In particular: 

w £ (x) = 1 + o(l) ondT £ r\D (14) 

Note that with this definition of /i^ fe , we have h £ ak (x) = e 2 h a ^ k (x / e) for 
n > 3 and h £ a k (x) = e 2 h ajk (x/e) + re 2 log e for n = 2. 

Proof. For every a, we denote by vf a the function 

vf a (x) = e 2 v a>£ -i D (x/£), 
defined in D and satisfying w £ a = on dD. 



if n > 3 



ao, ,,2 | e n r(k,u;) n 2 
2n p — eK| n_/ 

— la; — ek\ 2 — r(k,u))e 2 log \x — ek\ if n = 2, 
2n 
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1. For every a > ao, we have 

A(w £ -w £ Q )>a -a 
and Wq — w £ a = on dD. This implies 

wI(xq) - w £ a (x ) < / G(x ,x)(a - a) dx 



D 

where G(-,-) is the Green function on D (AG = 8 Xo and G = on 
dD). Note that we have 

G(xq,x) > —h(x — xo) Vx, xq G D, 

and so 

Wq(x ) - w £ a (x ) < (a - a ) / h(x-x )dx. 

Jd 

We deduce 

S uv(w £ -w £ )<{ C\D\y^ PD \a-a \ if n > 3 
supK u, a )<| cp | pDlogpi3 | Q _ ao | ifn = 2j 

with 

= inf {p ; L> C B p }. 

Hence we have 

Wq < w £ a + 0(a - q ). 
Using Lemma 03] (ii) (since a > ao), we deduce: 

w £ < h £ a k (x) + 0(a - a ) + o(l) Vx G B e/2 (ek) a.e. w G 0. 
which gives the second inequality in (p~3|) . 
2. Similarly, we observe that for every a < ao, we have 
A(w £ a -w £ )>a-a - al {tu , =0} . 
Proceeding as before, we deduce that for n > 3, 
sup « - i^) < Cp D llD^^iao -a) + Ca\{w £ a = 0}^^ 

D L 

and a similar inequality for n = 2. Using Lemma 14.31 (i), we get 

> K,k ~ o(e 2 ) - O(a -a)- Ca\{w £ a = O}^^. 
Finally, since 

lim|{< = 0}|=0 
for all a < ao, and (fl"3j) follows. 

□ 
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4.3 Proof of Proposition 12.21 

We are now in position to complete the proof of Proposition 12.21 We define 

w > 1 on T £ n D 



w e (x, u) = inf iw(x) ; Aw < a in D \ T £ , w = Q on dD \ Te 
it is readily seen that 

w £ (x,uj) = 1 on dT £ , 
< Aw £ (x, to) = ao on D\T e , 
. w s (x,uj) =0 on dD\T £ . 

So in order to complete the proof, we only have to show that w £ — > 
in H 1 (D)-we&k as e goes to zero. More precisely, we will show that w £ 
converges to zero in L p strong and is bounded in H . 

Strong convergence in LP: 

First of all, JH]) yields 



Wq{x) — o(l) < w £ (x,uj) < Wq(x) + o(l) ViGD £ a.e. well, 
which in turns imply (using Lemma 14.71 again) : 

K>k{ x )-°i l ) < w £ {x,lo) < h £ ak (x) + o(l) Vx G B £ / 2 {ek) a.e. u> € $7. 

(15) 

Next, a simple computation shows that 



f Ce n (e~+e 2 P) if n > 3, 

/ K, k \ p dx<{ v ; 

JB s \B a e y Ce 2 e 2v {\oge) p if n = 2 

Since #{eZ n n D} < Ce n for all n, we deduce from (fT5j) that 

f C fe^T +e 2N ) if n > 3, 
||w £ || L p < < V J ~ (16) 

[ Ce 2 (loge) ifn = 2. 

In particular 

w £ — ► in L p — strong, for all p G [1, oo). 
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Bound in H 1 : 

First of all, a simple integration by parts together with the fact that w £ = 1 
on dT £ yields 

/ \Vw £ \ 2 dx <a \D\+ [ \Vw £ \da(x) 

JD e JdT e 

where dT £ = UdS £ (k,uj). So we need an estimate in Vu> £ along dS £ (k,uj) = 

<9-B a =(r(fc,w))- 

We consider the function 



z(x) = < 



' w £ (x) - h £ a k (x) + |4 r 2 e n/(n-2) w h en n > 3 



It satisfies 



w £ (x) - h £ a k (x) + T^r 2 e 2 r when n = 2 



Az = in B 1/2 (ek) \ J B ae(r(fc ^ )) (e/c), 

z(x) = o(l) in B 1/2 (ek) \ B^^^ek) 
,z(x)=0 along dB a e {r{k;Ul)) (ek), 



and so 



n(n— 1J 

o(r n - 2 £ n e"^H = o(£ n a £ (r)^ n -^) if n > 3, 

if n = 2. 



|V*(s)| < 



(eV ^ 2 ) = o(e n a £ (r)- (n -^) 
on dB a e( r ( kuJ ))(£k). It follows that 

|V^| < |V/£, fc (s)| + |Vz(x)| < Ce n a e (r(k,u))-( n -V 

along dB a e( r ( kjU) )){ek), 
We deduce 

/ |Vu> £ | 2 dx < a |D|+ f \Vw £ \da(x) 

JD £ JdT e 

< a \D\+ / \Vw £ \da(x) 

< ao\D\ + Ce- n a £ (r) n - 1 e n a £ (r)'^ 

< c, 



and the proof is complete. □ 



20 



a £ (r) = re n/{n - 2 \ r(k,u) = ( - 

V 1 



5 Proof of Proposition 12.2b General case 

In this section, we treat the case where the sets S E (k,uj) have unspecified 
shape, but satisfy Assumption 2: 

cap{S e (k,u)) = e n j{k,uj). 

Throughout this section we assume n > 3. 

The proof makes use of the result of the previous section, after notic- 
ing that away from ek, the hole S £ {k,ui) is equivalent to a ball of radius 
a £ (r(k, uj)), where 

7(fc,*) y /(w - 2) 

n(n - 2)u> n J 

More precisely, we will rely on the following lemma: 

Lemma 5.1 For any k € Z n and uj € £1, let <p%(x,uj) be defined by 

f f v(x) > 1, Vie S £ (k,uj) ) 

<p%( x >v) = inf { v ( x ) ; Aw ^ > ,. /n „ r 

[ [ hm^i^oo = J 

T/ien /or are?/ 5 > 0, there exists R$ such that 

\(pl(x,u) - e n ^{k,ui)h{x - ek)\ < 5e n h(x - ek) 

for all x such that \x — ek\ > a e (Rs) and for all e > 0. 

Moreover, R$ depends only on the constant M appearing in Assump- 
tion 1. In particular, R$ is independent on k and uj. 

1. For a given 5 > 0, Lemma [5.11 implies that for every k € Z n and uj € f2 
there exists a constant R$(k,uj) such that 



(17) 

for all e > 0. Moreover, it is readily seen that for any R there exists 
£\{R) such that 

a £ (R)<e a /A for all s < e x . (18) 
for some a > 1. Finally, we note that by definition of cpt, we have 

|V<4| 2 dx = cap(£ e (fc)) = e n 7 (&,u/) (19) 
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2. Next, let ao and w £ be the coefficient and corresponding corrector 
constructed in the previous section, and associated with holes S £ of 
radius r(k,uj). Lemma |4. 71 implies that for 5 and R given, there exists 
62(6, R) < Si(R) such that for all e < 82(6, R), we have 

in B E/2 (ek), (20) 

in dimension n > 3. Note that thanks to (|18|) . Inequality (|2U|) holds 
in particular in B 2a ^ R ) \ B a e {R ){ek). 

The corrector given by Proposition 12.21 will be constructed by gluing 
together the functions <p e k (near the holes S e (k) and the function w £ (away 
from the holes). The gluing will have to be done in a very careful way so 
that the corrector satisfies all the properties listed in Proposition 12.21 For a 
given s, we define 5 £ to be the smallest positive number such that fll8|) and 
(|20p hold with 6 = S £ and R = R$ e . From the remarks above, we see that 
5 £ is well defined as soon as e is small enough (say smaller than £2(1, i?i)). 
Moreover, for any 5 > 0, there exists £0 = £2(81 Rs) such that 

5 £ < 5 V£ < £ . 

In particular 

lim 5 F = 0. 

From now on, we write 

Re = Rs e ■ 



We are now ready to define the corrector w £ : Let rj E (x) be a function 
defined on D such that 

rj e (x) = 1 on D \ U i£ z»B 2 a s (fl £ )(^) 
7] £ (x) = on U ke3d n B a e {RE) {ek). 

and satisfying 

|Vt? £ | < C^iReY 1 and \Ar) e \ < Ca £ {R £ )- 2 
in B2a^ Re ) \ B a E( Re }(k). We then define w £ (x,uj) in D by: 

w £ {x,uj) = r/ £ (x)w £ (x,uj) + (1 - rj £ (x)) ^ <p%(x, u) l Be/2 (ek)(%)- 
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w £ {x) 



e n r(k,uj) n - 2 
\x — ek\ n ~ 2 



It satisfies 

f V %(x) mB 2ae(Re) (k)\S £ (k) Vk£Z n 

W {X, LO) = < 

[ w £ (x) inD\U keZ »B a c {Rs y 
To simplify the notations in the sequel, we denote 

ip £ ( X ):= { Pk{ x ^) l B e/2 {ek){x) 

The properties of w £ are summarize in the following lemma, which im- 
plies Proposition 12.21 with (J3|) instead of the first equation: 

Lemma 5.2 The function w £ satisfies 

(i) uf = 1 on S £ for any e > 0. 

(ii) w £ converges to zero as e goes to zero in L P (D) strong for all p € [2, oo) 
and 

2n 

\\w E \\ LP < CepC»- 2 ) Vp > 2 
(Hi) w £ is bounded in H l (D). 

(iv) Aw £ converges to ao in ^(D) and thus satisfies 
Proof: 

(i) Immediate consequence of the definition of uf since </?| = 1 on S E (k, to). 

(ii) Assumption 1 yields 

(p e k (x,uj) < Ce n ~f(k,uj)h(x - ek) 

for all x such that \x — ek\ > a £ (M). Since tp £ k < 1 in B a er M \(ek), we 
deduce: 



v 

LP(UB R(k)a(s) ) 



(l-%)^ e Hi P(R n) < E fk 1 B E/2 (ek) 

k& n r\e- x D 

< E / 

+C E / (e n 7(fc)/i(x 

, _1 



< e « £ ( M ) r 



+C7 E e p "(« £ (« £ ))" _p(n " 2) 

k^Z n ne- 1 D 
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Using (fTBj) and the definition of a(e), we deduce: 



l-Ve)f e \\ P LP(Kn) < Ce- n M n e^+C^ Y e pn e n " p(n - 2) 



k<=z n nD 

fcez™nD 
< CM n e^ + C^e 2p 



where 2j? > if p > 2 and n > 3. 
Using (|16p . it follows that 



/ 2n \ 1/p 

IKHlp(d) < \\w \\lp(d) + C \£ n - 2 ) 

2n 

< Cep("- 2 ) 

for all p > 2. 

(hi) Next, we want to show that w £ is bounded in H 1 (D £ ). First, we note 
that in B £ / 2 (£k), we have: 

Vw £ = Vrj £ (w £ - tp%) + r, £ Vuf + (1 - r/ e )V<4 (21) 

where the function Vrj £ is supported in B 2a e^ R ^{ek) \ B a ^ R ^(ek) and 
satisfies 

iVTfcl <C(a e (R))-\ 
Since \w £ - <p%\ < in B 2a e^(ek) \ B a ^ Re) (ek), we deduce 

\V Ve {w £ -vi)?dx 

)(ek) 



f \V V£ (w £ - p £ )\ 2 dx < Y [ 



< Y (a £ (Re)T(a £ (R e ))- 2 ): 



k€eZ n nD 

< Ce- n e n = C, 



2(n-2) 



since we can always assume that 5 £ < 1 and R £ > 1. Finally, since w £ 
and ip £ are both bounded in H 1 (thanks to (|19p ). (|2ip implies 

||Vu; e || L 2 < C. 
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) It remains to evaluate the Laplacian of w £ . We have: 

Aw £ = a - (1 - Tj e )a + 2Vr] £ ■ V(w £ - (p £ ) + Ar/ e (w £ - cp £ ) in D E . 
Moreover, (J2DJ) and JT^) yield 

tie 

and by definition of w £ and tpf , we have 

A ( w£ ~ Lp k~^ X ~ efc | 2 ) = in B ^(R) \ B a*(R s )/2- 

Interior gradient estimates thus implies 

tie 

m B 2a-(R £ ) \ B a e (Re)- We deduce (using (USD): 

/ | Auf — a\ dx 
J D e 

< / (l-r) 6 )adx+ / |V%||V(w e - ip £ )\dx 
+ / |A%||t(; e — (p £ \ dx 

< E a£ (^) n 

+ £ ^(Rc)- 1 [ \V(w £ -vt)\dx 

+ e ae ( R er 2 [ \<-<p%\dx 



<c e « £ (^r+c E ^ aE (^rV(# 

< Ce {a ~ 1)n + C5 £ . 
In particular, 

lim / \Avf — ao\dx < C5 
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A Proof of Lemma 15.11 

We recall that n > 3 in this section. For any k G Z n , we define S £ (k) = 

_ n 

e n - 2 S e (k). Then Assumption 2 yields: 

cap(S e (£;)) = 7 (fe) <j. 

and Assumption 1 gives 

S £ (k) C B M (k). (22) 
For the sake of simplicity, we take k = 0. We recall that h is defined by 

h(x) ' ' 



n(n — 2)u) n \x\ n ~ 2 
Lemma 15. II will be a consequence of the following lemma: 
Lemma A.l Let (p be defined by 



ip{x) = inf < v (x) ; Av < , 



v(x)>l, Vx G S e (k, lo) 
\m\\ x \^ 00 v(x) = 



Then for any 5 > 0, there exists R, depending only on 5 and M such that 

\<f(x, oj) — jh(x) | < 5h(x) 

for all x such that \x\ > R. 

Proof: We recall that ip solves 

Aip(x) = for all x G W 1 \ S 

(f(x) = 1 for all x G S 

l im |x'Hoo fix) = 0. 

In particular, (|22p and the maximum principle imply 

M n ~ 2 

ip{x) < M n ~ 2 n{n - 2)oj n h(x) = in R n \ B M (0). (23) 

Next, we observe that 

= — / tpAtpdx= / \Vip\ 2 dx — / ifip^da(x) 
Ju n \s Jm. n \s JdS 
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and so 

|2 



Vip\ z dx = / fip u da(x) = / f u da(x). 
n \S JdS JdS 

Moreover, for any R > M, we have 

= / Atpdx = / c/?„ da(x) + / ip u da(x). 
JB R \S JdS JdBf, 



We deduce: 

7 = 



/ \V<p\ 2 dx = - <p u da{x) for all i? > M. (24) 



We now introduce the function 

6(x) = /i (h0 99 fe) =ri(n " 2K WF^ c 'fe) 

defined for x £ B 1 / M (0). A straightforward computation yields 

A9 = in B 1/M (0) 

and (J2"5j) implies 

G(x) < M n ~ 2 n{n - 2)u n in S 1/Af (0). 

A more delicate computation, making use of the mean formula for harmonic 
functions, gives 

f ip u da(x) = -Q{0). 

JdB R 

Hence fl2U) yields 

6(0) = eap(S e ) = 7 

To conclude, we note that interior gradient estimates for harmonic func- 
tions imply the existence of a universal C (depending only on M) such that 

|0(sc) -7| < C\x\ for all \x\ < 1/(2M). 

Inverting back, we deduce 

C 

\cp(x) - jh(x)\ < T-h(x) for all Id > 2M, 
\x\ 

which yields the result. □ 
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